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INTRODUCTION 

i 

kJj The action of a generalized fourth-order gravity is given by (for more details, see 01 and 

references therein): 

i* = \ J d 4 x^g [f(R) + 2^ m ] , (1) 

where ££ m represents the matter contribution, and the generalized field equations read 

00 1 ~ m R 

m ; G ab = T ab + T ah = T ab , (2) 

o 
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where 

Till 1 

T a m b = JTi T£ = j[W~ R f')Sab + Vftaf ~ gabVXf] • (3) 

Here we have defined /' = df(R) /dR and T% = \i m u a u b + p m h ab + q™u b + q^u a + 7t™ b , 
^ ; where ji m , p m , q m and K™ b denote the standard matter density, pressure, heat flux and 

anisotropic stress respectively. 

The total thermodynamics of the matter-curvature composite is then given by 

,,m yjn n m jj-m 

M-^ + M*, ^7+/> la^f+fi* = + (4) 

where [i R , etc. are thermodynamical quantities of the curvature fluid defined in the next 
section. The covariant derivative of a timelike vector u a can be decomposed into basic 
parts as 

V a u b = -A a u b + \h ab ® + o ab + £ abc co c , (5) 

where A a = u a is the acceleration, = V a w a is the expansion, o ab = V/ a u b \ is the shear 
tensor and co a = e abc V b u c is the vorticity vector. For the Weyl curvature tensor one has 

E a b = C abc dU C U d = E^afy , H ab = jEacdCfcU 6 = H( a fy , (6) 



giving a covariant description of tidal forces and gravitational radiation respectively. 



LINEARIZED FIELD EQUATIONS 

We consider the background to be Friedmann-Lemaitre-Robertson-Walker (FLRW), 
where the Hubble scale sets the characteristic scale of the perturbations. In the per- 
turbed spacetime the standard matter is considered to be a perfect fluid with the energy 
momentum tensor given by: 

T2 = {v m + P m )u a u b + p m g ab . (7) 

with p m = wfi m and the heat flux (q%) and the anisotropic stress (it™ b ) vanishing in the 
perturbed spacetime. In addition, since we consider shear-free perturbations, the shear 
tensor o ab vanishes identically. 

For the curvature fluid the linearized thermodynamic quantities are given by 

V R = j[h(Rf'-f)-®f''R + f"v 2 R] , 

P R = j, [\{f-Rf) +f''R + f"'R 1 + l{®f"R-f"y 1 R)] , 

<£=-Ji [f'RVaR + f'VaR - IfVVaR] , = J,f"V {a V b) R . (8) 

With the conditions above, the propagation and constraint equations can be given by 



0-V fl A fl = -l0 2 - ±(ju + 3p), (9) 

(co {a) )- - ±£ abc VbA c = -§0O) fl , (10) 

E (ab) _ e cd(ay cH )b = _ Q£ ab _ l^ab _ ^ q b ) _ 1 Qjtf ^ (1 1} 

H (ab) + eCd(ay cE )b = _ &R ab + \_ £ cd{ay ^b ^ ^ (U) 

fl m = -(flm+Pm)®, (13) 

A + V^f = -(n+p)®; (14) 

{Co) ab :=E ab -^ a A b ^-\nf = 0, (15) 

(d) fl := V fl - \e abc V b co c - \tf R = , (16) 

(C 2 ):=V a co a = 0, (17) 

(C 3 ) ab :=H ab + V^ a Cd^ =0, (18) 

(C 4 ) a := V a Pm + (n m + Pm )A a = , (19) 

(C 5 ) a := V h E ab + ±Vrf - \V a v + \®q a R = , (20) 



(C 6 ) fl := V b H ab + (n+ p )cQ a + \e abc V b q R c = 



(21) 



The conditions o ab = and q a m = give the two new constraints (Co) ab and {Cx) a 
respectively. Substituting (Co) b d into (Cs)^ and using (C4)/, we obtain the constraint 



w 



w + 



(22) 



where = lnjU m . To check the spatial consistency of the above constraint on any initial 
hypersurface we take the curl of (|22)) to obtain 



(O 



3/' . 



+ 



fw© I V 2 £o a = 



(23) 



where R = 2 (/1 — ^O 2 ) . Now defining the expansion, acceleration, jerk and snap pa- 
rameters by the following relations 



^ a da a 2 d 3 a 

= 3-, q = — ^, j 



a 3 dt 3 



a 3 d A a 
¥~d¥ 



(24) 



and using 



® = -^® 2 (l+q), q = - l -<d{j-q-2q 2 ) , = ^0 3 (2 + 3^ + j) 



j = l®( s + 2 j + 3qj) , q 



1 



-® 2 [s + 2j-3q 2 + 6qj-6q 3 ] , /? = -0Q 



(25) 



where 



£ = t© 2 (./-4-2)+£, 2 = i0[(4 + 5^ + j + ^ + 5)0 2 + 6^] 



(26) 



we can rewrite (l23l) as 



-+^-e 1 am 

2 3/'^' 



+ 



^2 + — 

./' 2. 



V 2 fti a 







(27) 



Spatial consistency requires the vanishing of either or the terms in the curly brackets. 
For temporal consistency differentiate ([27]) w.r.t. time to get 



0(O fl 
where 



\-w)P ~ (l+w)(3w + 5)/ + 4/"Q „ 

K H T, ~ Mm 



6/ 



Z 
+ P 



1 f fill fll \ fll 
Z =3[ l f-ty 2 ) 2 2 + ^7((4 + 5^ + j + ^ + ,)0 2 + 6 J R). 



0, (28) 



(29) 



It follows that for the new constraints to be spatially and temporally consistent we must 
have either ®oo a = or the expression in the curly brackets must vanish. It is easy to see 



from (|28l) that if the 3-curvature vanishes, then &co a = for vacuum universes (/i m = 0). 
This implies that a shear-free, spatially flat vacuum universe in any f(R) theory can 
rotate and expand simultaneously in the linearized regime. 

In the non - vacuum case, there exists at least one non-trivial case which does violate 
the Ellis condition. For a flat Milne universe, the Friedmann equation is given by 

2 d*f(R) f(R) Mo _ Q f3m 

and has the following general solution: 

„ . l+y^ l-y/3 Alln 3(l+w) 



Considering the particular solution {C\ = = C2), and comparing it with Eqn. (1281) . for 
the corresponding flat Milne universe in R" gravity, we obtain 

(1+Vv)Mm [3w + 9-4n]=0. (32) 



6/' 

Comparing solutions (1321) and the particular solution of (I3TI) (with n = 3(1 +w)/2) we 
find that w = 1 if \l m 7^ 0. In other words, for a stiff fluid in R 3 gravity, there exists a 
flat Milne-universe solution which can rotate and expand simultaneously at the level of 
linearized perturbation theory. 



DISCUSSION AND CONCLUSION 

In this work we showed that if the 3-curvature vanishes, then the result of |Q] can always 
be avoided for vacuum universes. We also demonstrated there is at least one physically 
realistic non-vacuum case in which both rotation and expansion are simultaneously 
possible. This suggests that there are situations where linearized fourth-order gravity 
shares properties with Newtonian theory not valid in General Relativity. 
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